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Partition density functional theory is a formally exact procedure for calculating molecular proper- 
ties from Kohn-Sham calculations on isolated fragments, interacting via a global partition potential 
that is a functional of the fragment densities. An example is given and consequences discussed. 



PACS numbers: 

Kohn-Sham density functional theory (KS-DFT) 0, Q 
is an efficient and usefully accurate electronic structure 
method, because it replaces the interacting Schrodinger 
equation with a set of single-particle orbital equations. 
Calculations with several hundred atoms are now routine, 
but there is always interest in much larger systems. Many 
such systems are treated by a lower-level method, such as 
molecular mechanics, but a fragment in which a chem- 
ical reaction occurs must still be treated quantum me- 
chanically. A plethora of such QM/MM approaches have 
been tried and tested, with varying degrees of success Q. 
These are often combined with attempts at orbital-free 
DFT, which avoids the KS equations, but at the cost of 
higher error and unreliability. 

On the other hand, partition theory (PT) 4, 5] com- 
bines the simplicity of functional minimization with a 
density optimization to define fragments (such as atoms) 
within molecules, overcoming limitations of earlier ap- 
proaches to reactivity theory [3, 0] • While there are now 
many definitions of, e.g., charges on atoms, none have 
the generality of PT and the associated promise of uni- 
fying disparate chemical concepts. However, previous 
work on PT has been either formal [3,0 or f° r t wo atom 
systems Hi]. 

In this paper, we unite KS-DFT with PT to produce 
an algorithm that allows a KS calculation for a molecule 
to be performed via a self-consistent loop over isolated 
fragments. Such a fragment calculation exactly repro- 
duces the result of a standard KS calculation of the entire 
molecule. We demonstrate its convergence on a 12-atom 
example. This also shows that fragments can be calcu- 
lated 'on the fly', as part of solving any KS molecular 
problem. 

Thus we present a formally exact framework within 
which existing practical approximations can be analyzed 
and, for smaller systems, compared with exact quantities. 
In practical terms, our method suggests new approxima- 
tions that can, by construction, scale linearly[lO| with the 



number of fragments (so-called O(N)), and allow embed- 
ding of KS calculations within cruder force-field calcula- 
tions (QM/MM). It also suggests ways to improve XC 
approximations so as to produce correct dissociation of 
molecules 

To understand the relation between DFT and PT, re- 
call that the Hohenberg-Kohn theorem proves that for a 
given electron-electron interaction and statistics, the ex- 
ternal (one-body) potential v(r) is a unique functional of 
the density n(r). The total energy can be written as: 

E[n]=F[n] + J d 3 rn(r)u(r), (1) 

where F[n] is a universal functional, defined by the Levy- 
Lieb constrained search [12] over all antisymmetric wave- 
functions W yielding density n(r): 

F[n] = min (*| f + V ee (2) 

*— >n(r) 

where T and V ee are the kinetic energy and Coulomb 
repulsion operators respectively. The KS equations are 
single-particle equations defined to reproduce n(r). De- 
fine the KS energy as 

E s [n] = ($ B [n]\f + V\<f> s [n])=T s [n] + J d 3 r n(r) v(r), 

(3) 

where V is the external potential operator, $ s [n] is the 
KS wavefunction (usually a single Slater determinant) of 
density n(r), and T s [n] is the kinetic energy of $s[n]. 
Define the Hartree-exchange-correlation energy, as 

E HXC [n] =E[n] - E s [n], (4) 

so that the KS potential is u(r) + v HXC (r), where 

«Hxc(r) = SE HXC [n]/Sn(r). (5) 



Partition theory deals with the problem of dividing a 
system into localized fragments. For molecules or solids, 



u ( r ) = E 
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v a (r), 



(6) 



where Zp is the atomic charge of a nucleus at point 
and these are regrouped into Nf fragment poten- 
tials, v a (r). The fragmentation is chosen based on the 
particular use of PT: e.g., one might atomize an entire 
molecule, or merely separate off a well-known chemical 
species. The partition problem is then to divide n(r) be- 
tween the fragments. There are many methods for doing 
so, but in PT[J,0] we minimize the total fragment energy 



E f 



min \^ ( F[n a ] + / d 3 rn a (r)v a (r) 
{"<*} \ J 



J2 a Li n Q (r)=n(r) 



(7) 

where n a (r) is the density on the a-th fragment. Each 
fragment is considered to be in contact with a distant 
reservoir of electrons [3], so its integral N a need not be 
an integer. If N a = p a + u a , with < v a < 1, then[l3T] 



where 



F[n a ] = (1 - v a )F[n Pa \ + v a F[n Pa+1 ], (8) 



n a (r) = (1 - is a )n Pa (r) + v a n Pa+1 (r), (9) 



and the integer densities are ground-states of a common 
potential. Minimizing the Lagrangian: 



G = E f + J d 3 r « p (r) Li(r) - £]n Q (r)j , 



(10) 



yields the correct fragment densities. The Lagrange mul- 
tiplier v p (r) is called the partition potential, and satisfies 



SF[r 



6n a (r) 



+ v a (r) + v p (r) = /i. 



So n a (r) is the ground-state density of N a electrons in 
effective fragment potential i>/ jQ (r) = v a (r) + v p (r). 

Thus PT replaces a molecule of interacting fragments 
with an effective system of non-interacting fragments, 
and v p (t) is the analog of v HXC (r) in KS theory. Analo- 
gous to Eq. we define the partition energy as 



E p [{n a }} =E[n] - E f [{n a }} 



(12) 



which, via Eq. ([7|), is a functional of {n a ('r)}a=i...N 
Functional differentiation yields: 



v p (r) = SE p [{n a }]/Sn a (r). 



(13) 



Analogous to KS-DFT, once £ p [{n Q }] is given, we have 
a closed set of equations that yield the molecular density 




FIG. 1: Solid line: The exact spin-unpolarized ground state 
of 12 electrons in the potential of Eq. (|17[l . Dashed lines: The 
fractionally occupied fragment densities. By symmetry, the 
other half of the density is simply the mirror image of that 
shown. 



and energy at self-consistency. Next, we use a superscript 
to denote quantities evaluated in the physical dissoci- 
ation limit, where all bond lengths between fragments 
have been taken to oo while keeping intrafragment dis- 
tances fixed. In this limit, the fragments do not interact 
and their densities do not overlap, and E ^ is the sum of 
the truly separated fragments, each with density n« (r). 
Defining the relaxation energy E re i = E^ — Ef, we write 



E n — E r l 



E r 



E d 



E-E 



(0) 



(14) 



where Edi S is the electronic contribution to the dissoci- 
ation energy. Thus E p < for any bound molecule (by 
construction), is expected to be much smaller than the 
total energy (on the scale of chemical bonding), and van- 
ishes as the fragments are pulled apart. 

In a KS fragment calculation, the KS potential for the 
a-th fragment is found from Eq. (I13p in KS quantities: 



(11) v s j,a[n a ,n a }(r) = v s [n a ](r)+(v(r) + w H xoN(r) - v s [n](r)) 



( 15 ) 

where w s [n](r) = —6T s [n]/Sn(r), and n(r) — n a (r) + 
n a (r). This is the central result of this paper, as it gives 
the fragment KS potential for a pair of trial densities, 
n Q (r) and n a (v), in terms of quantities from KS-DFT. 

In Eq. P^)l . w s [n a ](r) is simply the KS fragment po- 
tential from the previous iteration, but u s [n](r) is the KS 
potential for a trial density for the whole molecule. Many 
methods exist for finding this[l4j]. We iterate (IBj: 



v y s (r) 



4"°(r)+7 n^(r) 



dm) ( 



(r) , (16) 



where n^" 1 * 1 (r) is the density found from potential vi" 1 ^ (r) , 
7 > is a constant, and riM^ (r) is the target density 
(sum of fragment densities from the ft'th PDFT iteration) 
whose KS potential we are trying to find. To find the 



I I I I I I I 

-18 -15 -12 

X 

FIG. 2: The exact partition potential (solid line) for the at- 
omized chain and the fragment potential for the last atom 
(dashed line). The ground state with an occupation of 0.77 
in this potential can be seen as the end fragment density in 
Fig 1. 

fragment occupations, note that at self-consistency, the 
chemical potentials of all the fragments will be equal. 
We choose M fe+1) = Ni k) - T (u^ -£ (fc) ), where T 
is another positive constant and p, is the average of the 
fragment chemical potentials, used in conjunction with 
Eq. © for the functional [II]. 

The starting point in PDFT is to solve the KS equa- 
tions for each isolated fragment, generating their self- 
consistent densities (r) and KS potentials, and a trial 
molecular density which is the sum of overlapping atomic 
densities. Then Eq. (fT6|) is iterated to find its KS po- 
tential, which completes the inputs for Eq. (fl~S"|) . and the 
cycle repeated to self-consistency, when E<n s is found by 
evaluating E p on the final fragments, and subtracting the 
relaxation energies via Eq. ([Fiji . 

To show that our algorithm converges, we performed 
a PDFT calculation on a 12-atom Id chain with 12 spin- 
unpolarized non-interacting fermions, with potential: 



o^l cosh 2 [a; + (a - 6.5)R] ' 

We chose complete atomization into 12 fragments, so 
we only ever solve one- or two-electron problems in 
a single well. Fig. [T] shows the atomic and molec- 
ular densities after convergence. The molecular den- 
sity is identical to that found by direct solution of the 
eigenvalue problem for the entire molecule, and dou- 
bly occupying the first 6 eigenstates, which are delo- 
calized over the entire molecule. We see a small alter- 
nation between higher and lower densities throughout 
the molecule. The fragment density occupations reflect 
this, being 0.77,1.13,0.98,1.06,1.02,1.04 moving inwards 
towards the center of the chain. In Fig. fJl we show both 
the partition potential and effective fragment potential 
for the last atom. The (not very large) v p (r) polarizes 
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FIG. 3: The convergence of the fragment occupation values, 
N a , during an exact PDFT calculation. 

the density toward the molecular center, and shifts the 
density inwards compared to a free atom. The partition 
potential continues throughout the whole chain, lowering 
each fragment potential in the bonding region between 
atoms. The depth of these troughs oscillates, reflecting 
the oscillation in occupations. In Fig. [3J we show the 
convergence of the occupation numbers to their final val- 
ues, after some initial oscillations. The total energy of 
the molecular system can be found via Eq. (|T2|) . We find 
E f = -5.888 and E p = -1.803 leading to E = -7.691, 
which is exactly that of the direct solution. Since Ep is 
—6, \E re i\ <C \E P \ <C \E\, as expected. 

Our calculation was in fact far more expensive than a 
regular KS calculation, because we invert the KS problem 
for each trial molecular density exactly. But the purpose 
here was not speed, but the calculation of exact parti- 
tion potentials for small molecules and simple solids. It 
produces the exact partition potential corresponding to 
a given KS calculation for the molecule. 

The many potential uses of PDFT are made clear by 
this example. In principle, Eq. (|15j) is exact, but requires 
the KS potential of the entire system and to deduce the 
energy at the end of the calculation, one needs 




(18) 

where AG[n a ] — G[J2 n a] — E G[n a ]- However, 
any local-type approximation makes the method O(N). 
Thus, all the attempts of orbital-free DFT, to find useful 
approximations to AT S \n] , have now a simple framework 
in which to be tested[l6|. Moreover, there are no formal 
difficulties arising from taking density variations within a 
fixed density, as the trial molecular density is simply the 
sum of the fragment densities, which are varied freely. Al- 
though the exact fragment T s and v s (r) would be known 
during a calculation, approximations for AT S would take 
full advantage of any cancellation of errors. For embed- 
ding calculations, a simple approximation would be to 



4 



treat the system plus some fraction of its environment 
(a border region) exactly, and all the rest approximately. 
Since the KS potential is typically near-sighted, such a 
scheme should converge rapidly. 

For the dissociation of molecules, one can also see how 
to ensure correct dissociation energies within PDFT: sim- 
ply constrain occupations to be those of the isolated frag- 
ments. For H2, we constrain the spin occupations on the 
fragments to be (0,1) and vice versa. Of course, this 
is what happens when symmetry is broken as the bond 
is stretched, and the difficulty is in producing a scheme 
that seamlessly goes over to (1/2,1/2) occupations as R 
reduces to the equilibrium value. The value of our formal- 
ism is that it produces a framework for both addressing 
these questions and constructing approximate solutions. 

There is a simple adiabatic connection formula for 
PDFT. Consider scaling all bond lengths between frag- 
ments by A -1 (again keeping intrafragment densities 
fixed), where < A < 1. For each A, we find those 
molecular densities whose fragment densities match those 
of our molecule, and define the corresponding partition 
energy, E p (X). At A = 1, we have the original molecule; 
as A — > 0, the bonds become large and the fragments 
do not interact, so that E p (0) — 0. For intermediate A, 
the molecular density is simply that of the fragments, 
overlapped a distance R/X apart. Then 



which produces bizarre functional derivatives, different 
from those of KS DFT. None of these issues arise once 
smooth (e.g. local or gradient-corrected) approximations 



E = Ef+ l 



dEp{X) 
dX ' 



(19) 



This allows all the methods of traditional intermolecular 
symmetry-adapted perturbation theory (SAPT) [13] to be 
applied to this problem, but with the advantage that the 
fragment densities remain fixed. Interestingly, because 
the fragments will generally have dipole moments, the 
partition energy decays as 1/i? 3 , so that the integrand 
above behaves as A 2 . (For physical systems that are well- 
separated and have attractive van der Waals forces, such 
effects must be cancelled by analogous terms in E re i). 

There has been considerable previous work on schemes 
designed to allow a fragment calculation of a larger 
molecule, either within the framework of orbital-free 
DFT or atomic deformation potentials, sometimes pro- 
ducing the same (or similar) equations. Among the ear- 
liest, Cortona's crystal potential (later called embed- 
ding potential) [H, EH is an intuitive prescription for 
v p (r). But our formalism reproduces the exact solution 
of the original problem, using only quantities that are 
already defined in KS-DFT. For example, this is not 
possible in general without the ensemble definition of 
Eq. dHJ), which produces the correct self-consistent oc- 
cupations (unlike, e.g., the self-consistent atomic defor- 
mation method [20l. |21| . where this choice leads to a ba- 
sis set dependence [22]) . We also never freeze the total 
density 23, 24, [25|, but only ever consider it as a sum 
of fragment densities. This avoids ever needing density 
variations that are limited by some frozen total density, 



are made to the kinetic energy functional[l8|, [19|, |21| . 
but they are vital in a formally exact theory. Thus the 
present PDFT can be regarded as a formal exactification 
(and therefore justification) of these pioneering works. 

KB and PE acknowledge support under NSF grant 
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